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We study the hidden sector effects to the mass renormahzation of a simplest gauge-mediated 
supersymmetry breaking scenario. We point out that possible hidden sector contributions render 
the soft scalar masses smaller, resulting in drastically different sparticle mass spectrum at low energy. 
In particular, in the 5 + 5 minimal gauge mediated supersymmetry breaking with high messenger 
scale (that is favored by the gravitino cold dark matter scenario), we show that stau can be the next 
lightest superparticle for moderate values of hidden sector self coupling. This provides a very simple 
theoretical model of long-lived charged next lightest superparticles, that imply distinctive signals in 
ongoing and upcoming collider experiments. 



I. INTRODUCTION 

The minimal supersymmetric standard model (MSSM) 
is an attractive framework of a theory beyond the stan- 
dard model. It resolves some of the shortcomings of 
the standard model and at the same time makes falsifi- 
able predictions for upcoming experiments. The presence 
of dark matter, for an instance, has been consolidated 
by numerous cosmological and astrophysical observations 
According to the best- fit A-cold dark matter (the 
'concordance') model, the major part of the energy den- 
sity of the present universe is dark energy and dark mat- 
ter, neither of them is explained by the standard modeL 
While dark energy is Ukely to be of gravitational origin, 
dark matter is readily accommodated in MSSM; as the 
Ughtest superparticle (LSP) is protected by R-parity, it 
is an ideal candidate of dark matter. 

Among competing supersymmetry breaking scenarios 
of MSSM, the gauge-mediated supersymmetry breaking 
(GMSB) is advantageous for natural suppression of 
flavor-changing neutral currents and CP violation. In 
GMSB the LSP is gravitino. Cosmological model build- 
ing with successful structure formation favors cold dark 
matter, which, in the case of gravitino dark matter, yields 
gravitino mass lower bound 771,3/2 ^ 100 keV. Gravitino 
LSP of this mass range is realized by GMSB models with 
sufficiently high messenger scale (see also Section 3). For 
this reason we will be concerned with the high messen- 
ger scale GMSB scenario in this paper. The interaction 
of gravitino with the standard model particles is Planck- 
suppressed. While this makes gravitino an ideal candi- 
date of dark matter, it also makes its direct detection 
rather difficult. A more direct telltale of supersymmetry 
is from the next lightest superparticle (NLSP), which. 
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in the case of GMSB, is either the lightest neutralino 
or the lightest charged slepton(s). In the high messen- 
ger scale GMSB, the NLSP is long-lived, behaving as 
quasi-stable particles, and the distinction of whether the 
NLSP is neutralino or slepton leads to entirely different 
physics. In the 'minimal' GMSB scenario based on the 
A^5 = 1 (the sum of the Dynkin indices) 5-1-5 messen- 
gers which are representations of the SU{5) gauge group 
[3, 01, the NLSP is neutralino. The long-lived neutralino 
NLSP scenario is strongly constrained by the big-bang 
nucleosynthesis (BBN) Q, leading to severe constraints 
on both gravitino and neutralino masses. Assuming non- 
thermal production of gravitino dark matter combined 
with the BBN bound, the neutralino mass is constrained 
to be larger than TeV scale. In contrast, if the gravitino 
dark matter is of thermal origin the neutralino mass can 
be smaller. In any case, the long-lived neutralino be- 
haves like stable particles in collider experiments and the 
physics is similar to that of neutralino LSP. If the NLSP 
is charged slepton(s), on the other hand, we expect en- 
tirely different signals. A natural candidate for a charged 
NLSP is the lightest scalar tau (stau). While the stau 
NLSP scenario is also constrained by BBN Q , the restric- 
tions on the gravitino and stau masses are not so tight 
as the neutralino NLSP case. In collider experiments, 
charged particles like stau leave tracks in detectors, al- 
lowing direct observations of the NLSP. It would also be 
possible to make a precise measurement of its mass, and 
the absence of missing energy (such as due to neutrino) 
would facilitate detailed study of various processes in- 
volving stau jel-fiot. In GMSB, stau NLSP is obtained 
in models with larger TVs 0] or messengers belonging 
to larger gauge groups [ll[. The messenger index 7V5, 
however, cannot be taken to be arbitrarily large as it is 
constrained by the condition that the successful gauge 
group unification must be preserved. In the case of sim- 
ple multiple SU{5) messengers, for example, the possible 
range of is roughly iVs < 10 for the messenger masses 
« 10^° GeV. For too large A^5 the gauge coupUngs diverge 
before the unification. 



2 



These studies are based on the traditional assumption 
that the renormalization group (RG) flow of the visible 
world masses does not depend on details of the hidden 
sector dynamics, which, as is realized recently [T2|. [l3j. 
turned out not necessarily to be the case. Based on this 
observation, effects of hidden sector on the renormaliza- 
tion group equations (RGE) are studied in [ij]- They use 
a toy model of the hidden sector with a self-interaction 
coupling in the constrained MSSM, and show that the 
soft mass spectrum at low energy is altered by hidden 
sector effects, suggesting that such spectra can be used 
to determine the hidden scale and the strength of the 
self-coupling. In the present paper we shall carry out 
a similar analysis on the GMSB scenario that gives rise 
to gravitino cold dark matter, by taking into account the 
contributions from the hidden sector and reexamining the 
RG analysis of the minimal GMSB d i]. We find that 
the effect of the renormalized self-interaction coupling 
of the hidden sector field renders the soft scalar masses 
smaller than the conventional GMSB results, and con- 
sequently the low energy mass spectrum of the MSSM 
particles can be drastically altered. A phenomenologi- 
cally interesting consequence of our analysis is that stau 
can be the NLSP even in the minimal GMSB scenario. 
Our model thus provides a simple stau NLSP model with 
gravitino cold dark matter. 

The plan of the paper is as follows. In the next sec- 
tion we present the setup of our model, focusing on how 
the hidden sector dynamics contributes to the RG of the 
visible sector. In Section 3 we present our numerical 
results, and in Section 4 we discuss constraints of the 
parameter space from the collider and cosmological ex- 
periments. We conclude in Section 5 with comments, and 
explicit forms of the hidden sector contributions to the 
RGE are summarized in the Appendix. 



II. THE HIDDEN SECTOR IN THE 
GAUGE-MEDIATED SUPERSYMMETRY 
BREAKING 

In the GMSB scenario, the messengers couple directly 
to the hidden sector fields (which is responsible for the 
supersymmetry breaking), and indirectly to the MSSM 
matter fields through the standard model gauge interac- 
tions. Once the supersymmetry is broken in the hidden 
sector, the gauginos and the scalars in the MSSM become 
massive via one-loop and two-loop corrections, respec- 
tively. In the minimal GMSB based on the A^s = 1 5 -f- 5 
messenger belonging to an SU{5), the gaugino masses at 
the messenger scale M are given by 



Mait = 0) 



where 
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t = ln(/i/M) and A = F/M with F the supersymmetry 
breaking scale. The index a = 1,2,3 is for the MSSM 
gauge groups C/(l), SU{2), SU{3), and ga are the corre- 
sponding gauge couplings. The soft scalar masses at the 
messenger scale are given by 



r, = ^ ^ 
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where the index i denotes the MSSM scalars 
{q,u,d,l, Hu, Hd) in this order, and C2{Ri) are 
the quadratic Casimir for the matter fields in represen- 
tation Ri of the a-th MSSM gauge group. As indicated, 
these masses are generated at the messenger scale M . 
Physical masses at the electroweak scale are evaluated 
through RG evolution. Below we follow [l^, [ij] and 
derive the RGE, taking the hidden sector effects into 
account. 

We model the hidden sector by a single chiral superfield 
X with Lagrangian, 



d^eW + h.c. 



where the superpotential is 



W^-X\ 
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This does not itself break the supersymmetry but we as- 
sume that the F-term component of X acquires a non- 
zero vacuum expectation value by some mechanism. We 
shall be interested in the effects of the hidden sector 
renormalization. The hidden sector field X couples only 
to the messenger fields, and the messenger fields couple 
to the MSSM fields through the gauge interactions. Be- 
low the messenger scale M, the messenger fields are in- 
tegrated out, yielding the following effective interactions 
between X and the MSSM fields. 



^XX^ 
~M2 



d^e—w^w^ + h.c. 



(6) 



with real and complex coefficients fc^ and Wa ■ Here (f>j are 
the MSSM matter superfields (their scalar components 
are (j), = q, u, d, I, e, Ha) and are the MSSM field- 
strengths. As the F-tcrm component of X obtains a non- 
zero vacuum expectation value at the messenger scale M, 
the first and the second terms of ([6]) respectively give the 
scalar masses ^ and the gaugino masses (P). 

Let us consider quantum corrections to There 
are no 1 particle irreducible diagrams that renormalize 
operators linear in X, and hence the second line of ^ 
receives only the wave function renormalization coming 
from the hidden sector field. By the renormalization 
X Zj^ X, the second line of ([6]) becomes 



^ M " 
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The integrand is a holomorphic function whereas the 
is real. Therefore, due to the non- 



factor Z 



-1/2 



X 



rcnormahzation theorem [isl there is no quantum 
correction in the second hue in ([S]). This means that Wa 
does not receive any quantum correction at ah scales. On 
the other hand, the integrand of the first term in ([6|) is 
not a holomorphic function and therefore the coefficients 
hi are rcnormalized by the visible sector gauge interac- 
tions and the hidden sector coupling A in ([5]). They then 
satisfy the RGE 



dt 
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The second term is the leading visible sector contribution. 
In the first term "/(t) is the anomalous dimension at t = 
ln(/.j/Af) arising from the hidden sector interactions ([Sj. 
The anomalous dimension in the lowest order in A is 



7W 
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where X(t) is the running hidden sector Yukawa coupling 
which renormalizes according to 



dX{t) 



The RGE (H 

ki(t) = cxp 
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are solved as 

•0 



dt'-f{t') h{Q) 
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Using this expression, the scalar masses of the visible 
sector including the RG effects are described by 

mf{t) = h{t)A^, 



with 
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Here the gauge couplings run according to the standard 
1-loop formula, 
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where (6f,6f,6f) = (-3,1,33/5) for the MSSM and 
(^^1,^2, 63) = (-7,-19/6,41/10) for the standard model 
gauge couplings. Mz and Mg are the Z-boson mass and 
a typical soft mass scale, respectively. 
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FIG. 1: The slepton mass RG flows with and without the 
hidden sector efi'ects. The upper curves (black and blue) are 
the left-handed, and the lower curves (red and orange) are 
the right-handed sleptons. The 1st and the 3rd generations 
are respectively indicated by the solid and the dashed curves. 
The curves with sharp decline (blue and orange) above the 
hidden scale Mhid = 10*'^ GeV (indicated by the vertical 
dotted line) are the flows including the hidden sector effects. 
The straight lines (black and red) are the flows without the 



hidden sector effects. We have chosen M 
and tan /3 — 10. 



10'^ GeV, A = 3.8 



III. THE HIDDEN SECTOR 
RENORMALIZATION GROUP FLOW 

We shall be interested in the low energy mass spec- 
tra, in particular determination of the NLSP. The hid- 
den sector renormalization affects the mass RG flows be- 
tween the messenger scale fi = M and the hidden scale 
/i = Mhid- Eq. © implies that the RGE for the MSSM 
soft scalar masses are modified as. 



"df 



dmf AAt 



dt 
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while those for the other masses and couplings are un- 
altered (see, for example, [13 )■ Below the hidden scale, 
the hidden sector fields are integrated out and all the 
masses and the couplings evolve according to the stan- 
dard MSSM RGE. The low energy mass spectrum is ob- 
tained by integrating the RGE, from the messenger scale 
to the hidden scale with the hidden sector effects 
included, and then down to the electroweak scale in the 
standard way. In this section we present numerical study 
of the MSSM RGE including the hidden sector effects at 
one-loop order. The non-standard parts of the RGE are 
listed in Appendix [X] 

In our computations we made following approxima- 
tions. For simplicity only the (3, 3) family component 
of the three Yukawa matrices, and only the (3, 3) family 
component of the trilinear A-term matrices, are set to be 
non-zero. The latter are set to vanish at the messenger 
scale. Throughout our analysis we use following values: 
tan/3 — 10, the soft mass scale Ms = 500 GeV, and 
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FIG. 2: The squark mass RG flows with and without the 
hidden sector effects. The curves starting from the larger 
mass at M = 10^^ GeV are the left-handed squarks, and those 
starting from the smaller mass are the right-handed up-type 
squarks. The 1st and the 3rd generations are denoted by 
the solid and the dashed curves. The blue (left-handed) and 
orange (right-handed) curves with sharp decline above the 
hidden scale Mhid = 10*'^ GeV (the dotted vertical line) are 
the flows with the hidden sector effects, whereas the black 
(left-handed) and red (right-handed) curves show the flows 
without the hidden sector effects. The parameters are the 
same as in Fig. [1] (A = 3.8 and tan /3 = 10). 
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FIG. 3: The stau/neutralino mass ratio mf/m^o plotted 

against the hidden sector coupling A. We used M = 10^^ 
GeV and tan/? = 10. 



A = FjM = 10^ GeV. The hidden sector scale is defined 
as Mhid = = VA X M. For instance, for M = 10^^ 
GeV, the hidden scale is Afhid = lO^'^ GeV. 

Fig. [T] shows the mass RG flows of the flrst and the 
third generation slcptons with and without the hidden 
sector effects, and Fig. [2] shows similar results for the 
squarks. The messenger scale is taken to be M ~ 10^^ 
GeV. In both Fig. [T] and Fig. [2l we see that the hidden 
sector renormalization (down to the hidden scale Miiij) 
gives rise to smaller soft scalar masses, compared to the 
flows without it. This effect is readily understood by 



looking at the RGE ((M|) - ((J7| . in which the effects of 
the hidden sector renormalization are the positive terms 
4A^m| on the right hand sides^. It is also obvious from 
the RGE that larger A gives smaller soft scalar masses 
at low energy. In contrast, the RGE and hence the low 
energy mass spectrum of the gauginos are not affected 
by the hidden sector renormalization. This implies that 
the lightest scalar supcrparticle, namely the lightest stau, 
can be lighter than the lightest gaugino, i.e. bino, for suf- 
flciently large A. In Fig. [31 we plot the stau/neutralino 
mass ratio mf/m-^o at the electroweak scale as a func- 
tion of A with M = 10^^ GeV. A = corresponds to no 
hidden sector RG effects, reproducing the known result 
of neutralino NLSP in the minimal GMSB scenario. As 
A is increased^, the mass rrif of the lightest stau f be- 
comes smaller, while the neutralino mass m^o remains 
the same; the stau becomes lighter than the neutralino 
when A > 3.7. For example, M = lO^^ geV and A = 3.8 
yields rrif = 132 GeV and the neutralino (bino) mass 
m^o = 135 GeV at the electroweak scale. This stau mass 
is compatible with the LEP bound for long-lived massive 
charged particles > 102.0 GeV [ll|. 

The ratio of the stau and the bino masses m^/m^o 
is depicted in Fig. 21 against the messenger scale M for 
various values of A. The vertical dashed line indicates 
the minimal messenger scale determined by the condition 
that gravitino is cold dark matter, m^/^ ^0.1 McV. Here 
the gravitino mass is given by [13, [2^ 



where Mp = 2.4 x 10^^ GeV is the reduced Planck mass. 
The gravitino mass is also affected by the hidden sector 
renormalization (TJ]. For example, when the source of 
the supersymmetry breaking is a term linear in X, the 
vacuum expectation value of the i^-term component of 
X is proportional to the wave function renormalization 
factor Z-^^^ = exp(— i dt^it)). This takes effect from 
the messenger scale down to the hidden scale, resulting in 

— 1/2 

suppression of the gravitino mass by the factor « 
0.5 for A = 2 ~ 4. Taking this suppression into account, 
the minimal messenger scale consistent with the cold dark 
matter condition mg/j > 0.1 MeV is k, 9.0 x 10^ GeV. 
By choosing an appropriate messenger scale M within 
< 10^2 GeV, stau can be the NLSP with A as 
small as 3.0. 



^ In GMSB a hidden sector field need be an MSSM gauge singlet 
and a possible model of hidden sector superpotential is Yukawa 
type, as l(5]l. Then the hidden sector contributions to RGE are 
always positive. 

^ The perturbative means is valid for A^/(47r) < 1. Our choice of 
the parameter in the following discussions is roughly within this 
range. 
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FIG. 4: The stau/neutralino mass ratio mf/m^o plotted 
against the messenger scale M, for A = 2 (solid red), A = 3 
(dashed brown), A = 4 (dot-dashed purple), A = 5 (dotted 
orange) from above, with tan /3 = 10. The vertical dashed 
line indicates the minimal messenger scale in our model, 
M, = 9.0 X l(f GeV. For large enough A stau becomes lighter 
than neutralino. 



IV. STAU NLSP PHENOMENOLOGY 

In the previous section we have seen that the hidden 
sector renormalization can change the low energy mass 
spectrum. Here we describe phenomenological implica- 
tions when the NLSP is the lightest stau. 

We note that cosmological constraints on the mass 
parameters in the stau NLSP scenario are less strin- 
gent than in the neutralino NLSP case Q. In BBN, 
stau produces a bound state with ''He that tends to 
overproduce primordial ^Li through catalytic reaction 
^He f- + D ^ '^Li f- HH, m (see also 0, H, (H), 
giving upper bound for the stau lifetime'^ 

4877771^ ,„Mp / mLoN ^ 
T, = T-\f ^ tG) 5^ 1 - ^ 

< 5 X 10^ seconds. (17) 

Fig. E] shows the stau lifetime Xf plotted against the stau 
mass TOf for different values of gravitino mass 7713/2- The 
upper bound of the stau lifetime Tf < 5 x 10^ con- 
strains the gravitino mass from above. For our typical 
range of stau mass 777^ 130 GcV, the gravitino mass 
is 7773/2 ^ 0.6 GeV. Together with the condition that 
gravitino behaves as cold dark matter, a successful cos- 
mological scenario with stau NLSP and gravitino cold 
dark matter requires the gravitino mass to be within 0.1 



FIG. 5: Lifetime of stau plotted against the stau mass, for 
= 1 MeV (red solid curve), 7713/2 = 10 MeV (brown 
dashed), 7713/2 = 0.1 GeV (purple dot-dashed) and 7773/2 = 1 
GeV (orange dotted). 

MeV < 7773/2 < 0.6 GeV. The model prediction is 
safely within this range. 

The origin of relic gravitino is cither thermal or non- 
thermal productions. The thermal production is due to 
scattering in the thermal plasma in the radiation domi- 
nated era after inflation. The relic abundance from this 
contribution is computed as (26l - [28j 

where Tr is the reheating temperature and h is the Hub- 
ble constant in the unit of 100 km Mpc~^s~^. The non- 
thermal production of gravitino is due to decay of the 
NLSP * [SSl [39} . The abundance of non-thermally pro- 
duced gravitino is estimated as (ioMisj 

G - ^ -"-^ VlOOGevAlTevJ' 

(19) 

where fi^' is the thermal abundance of stau. The total 
gravitino abundance is thus 

n^h^ ^n^h^ + nl^^h^ (20) 

which is constrained by the observed dark matter density 

a 

f^G^^ < f^CDA/^^ ^ 0.1131 ±0.0034. (21) 

For stau mass ?77f w 130 GeV, the BBN bound for 
the gravitino mass 7773/2 < 0.6 GeV restricts the non- 
thermally produced gravitino abundance ()19|) to be ex- 
tremely small, il.^^h'^ < 10^^, indicating that most of 



^ This constraint may be relaxed by considering substantial left- 
right mixing of the stau mass eigenstates [25II . 



* The productions of gravitino from decay of moduli, inflation, 
and Polonyi fields have also been discussed [29l437ll . 
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the gravitino dark matter is produced thermally. The 
BBN bound for the gravitino mass and the dark matter 
abundance ((2T|l gives an upper bound of the reheating 
temperature. Using we find Tr < 10^ GeV. The 

model prediction ([16]) with M = lO^^ Q^y and A ^ 10^ 
GeV yields 7713/2 « 11 MeV; if the right amount of grav- 
itino dark matter of this mass is to be produced thermally 
{flcDM = i^^^), one needs Tr w 10^ GeV. This is lower 
than the reheating temperature required for successful 
leptogenesis; the observed baryon asymmetry needs to 
be generated by some other mechanism, such as [44 1. 

Stau NLSP is also of interest for collider phenomenol- 
ogy. Eq. (17) implies that the lifetime of stau Tf may 
be long, depending on the gravitino and stau masses; for 
instance, Tf « 100 sec for nif = 130 GeV and 1113/2 = 0.1 
GeV. For a particle with such a long lifetime, the decay 
length well exceeds the size of detectors and the decay 
takes place outside detectors. There have been intrigu- 
ing proposals for trappingsuch long-lived charged NL- 
SPs outside detectors |6l-[ig|. It has been pointed out 
that detailed studies of stau decay may provide precise 
measurements of the gravitino mass and the supersym- 
metry breaking scale. This would also be a crucial test 
of supergravity since the stau decay f — tG involves 
supergravity effects. 



V. DISCUSSION 

In this paper we studied the effects of hidden sector 
renormalization in the minimal GMSB model, using a 
toy model of the hidden sector field X with self-coupling 
A. We saw that the soft scalar and gravitino masses are 
affected by the hidden sector RG; they become smaller 
than the conventional GMSB RG results. Since the gaug- 
ino masses are unaltered, the resulting low energy mass 
spectrum of the supcrparticles can be quite different from 
the conventional GMSB. In particular, we find that stau, 
instead of neutralino, can be the NLSP in the minimal 
GMSB model. Based on this observation we also dis- 
cussed phenomenological implications of stau NLSP and 
gravitino cold dark matter arising from the model. With 
our typical choice of parameter values: the messenger 
scale M = lO^^ GeV, tan ^ = 10, A = F/M = 10^ GeV 
and the hidden sector coupling A = 3.8, the model yields 
stau NLSP with mass rrif ~ 130 GeV. The gravitino 
mass arising from these parameter values is shown to be 
consistent with the BBN constraints. It was also shown 
that gravitino cold dark matter with observationally con- 



sistent abundance is attributed to thermal production in 
this scenario. 

We comment that the hidden sector dependence of the 
mass spectra, as discussed in this paper, does not nec- 
essarily mar the predictability of the GMSB scenario. 
Rather, the hidden sector only influences the soft scalar 
and the gravitino masses, and the beauty of the GMSB 
scenario is kept intact. In fact, as emphasized in [T2| - [T^ . 
the effects of the hidden sector RG can be viewed as a 
new window to look into the hidden sector using the low 
energy mass spectrum that is expected to be revealed in 
experiments such as Large Hadron Collider. While the 
model discussed here represents a simple toy example of 
many possible hidden sector models, we believe essential 
features of our results to be generic. For example, only 
the soft scalar and the gravitino masses, and not the other 
masses neither the couplings, are subject to the hidden 
sector RG effects. This is a consequence of the powerful 
non-renormalization theorem and we do not expect this 
feature to be dependent on details of the hidden sector. 
Needless to say, however, there might be various other 
mo del- dependent issues that deserve further study. 
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Appendix A: The renormalization group equations 

In this appendix we list the RGE for the soft mass 
parameters that involve the hidden sector contributions. 
The new ingredient is the hidden- visible interaction ([5|), 
which give additional terms to the soft mass RGE. 
Within our approximation only the (3,3) family compo- 
nent of the three Yukawa matrices, yt, ?/&, j/t, are set to be 
non-zero. The corresponding non-zero components of the 
three trilinear A-term matrices are at , , a,- . We use nor- 
malized variables At = at/uti ^6 = o-hlvbi = o-tIVt- 

The RGE for the soft masses are 



8^^ 



~di 
,dml 
~dF 



26 
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15 



(Al) 
(A2) 
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where 

6 = ytimH^+ml + ml + A't), (A8) 
6 - (™L +ml + ml+Al), (A9) 

and 

+ Tr(?7iQ - 2m^ + m|, + m| - ?7i|)} . (All) 



(A3) 
(A4) 
(A5) 
(A6) 
(A7) 



Here the trace means the sum over the generations. The 
above equations for the squarks and sleptons apply to the 
3rd family; the equations for the 1st and the 2nd families 
are obtained from above by simply omitting the Yukawa 
and the A-term contributions (i.e. neglecting the ^t, £,b, 
6 terms). All the other RGE (the gauge couplings, the 
gaugino masses, the Yukawa couplings and the A-terms) 
are unaffected and are given e.g. in |17| . 
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